arXiv: 1506.02704vl [astro-ph.SR] 8 Jun 2015 


Astronomy & Astrophysics manuscript no. v6 ©ESO 2015 

June 10, 2015 


Seismic evidence for a weak radial differential rotation in 
intermediate-mass core helium burning stars 

S. Deheuvels^’J. Ballot^’P. G. Beck^, B. Mosser"^, R. 0stensen®, R. A. Garcia^, and M. J. Goupil^ 

^ Universite de Toulouse; UPS-OMP; IRAP; Toulouse, France 
^ CNRS; IRAP; 14, avenue Edouard Belin, F-31400 Toulouse, France 

® Laboratoire AIM Paris-Saclay, CEA/DSM-CNRS-Universite Paris Diderot, IRFU/SAp, Centre de Saclay, 91191 
Gif-sur-Yvette cedex, France 

^ LESIA, Observatoire de Paris, PSL Research University, CNRS, Universite Pierre et Marie Curie, Universite Denis 
Diderot, 92195 Meudon cedex, France 
® Instituut voor Sterrenkunde, KU Leuven, 3001 Leuven, Belgium 

June 10, 2015 


ABSTRACT 

Context. The detection of mixed modes that are split by rotation in Kepler red giants has made it possible to probe the 
internal rotation profiles of these stars, which brings new constraints on the transport of angular momentum in stars. 
Mosser et al. (2012) have measured the rotation rates in the central regions of intermediate-mass core helium burning 
stars (secondary clump stars). 

Aims. Our aim is to exploit the rotational splittings of mixed modes to estimate the amount of radial differential 
rotation in the interior of secondary clump stars using Kepler data, in order to place constraints on angular momentum 
transport in intermediate-mass stars. 

Methods. We select a subsample of Kepler secondary clump stars with mixed modes that are clearly rotationally split. 
By applying a thorough statistical analysis, we show that the splittings of both gravity-dominated modes (trapped in 
central regions) and p-dominated modes (trapped in the envelope) can be measured. We then use these splittings to 
estimate the amount of differential rotation by using inversion techniques and by applying a simplified approach based 
on asymptotic theory (Goupil et al. 2013). 

Results. We obtain evidence for a weak radial differential rotation for six of the seven targets that were selected, with 
the central regions rotating 1.8±0.3 to 3.2± 1.0 times faster than the envelope. The last target is found to be consistent 
with a solid-body rotation. 

Conclusions. This demonstrates that an efficient redistribution of angular momentum occurs after the end of the main 
sequence in the interior of intermediate-mass stars, either during the short-lived subgiant phase, or once He-burning 
has started in the core. In either case, this should bring constraints on the angular momentum transport mechanisms 
that are at work. 
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1. Introduction 


Even though rotation is known to have an important im¬ 
pact on stellar structure and evolution, we still lack a the¬ 
oretical understanding of how the internal rotation rates of 
stars evolve in time. Several physical processes can trans¬ 
port angular momentum (AM) in star s: hydrody namical 
instabilitiesa£d_meridionalmirculatipnj[Zahn|ll99^Mathis 
& Zahn l2004l ). internal gra vity waves (IGW) excited a t the 
edge of convective regions (ICharbonnel fc Talonll200^ , and 
magn etic mechanisms, either fr om a buried fossil magnetic 
field (GouEfh fc Mclntvrd 19981) or from m agnetic instabil¬ 
ities ( Spruitl Il999l . iRudiger et ^ l2014lJ ). It is now well 
known that rotationally-induced mechanisms of AM as they 
are currently understood are far too inefficient to produce a 
rigid rotation for the solar radiative interior as found by he- 
liosei smology (ISchou et al.]ll998l.lChaplin et ffi]lI999l. Garcia 
et ai. l2nnsfi . 
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The asteroseismic measurement of the internal rota¬ 
tion profiles of red giants has recently shed some new light 
on this long-standing problem. The very high precision 
reached by the Kepler space mission has made it possible 
to measure the rotational splittings of mixed modes, i.e. of 
modes that behave both as pressure modes in the envelope 
and as gravity modes in the core. It was thus shown that 
stellar£ores_S 2 iiyi 2 i.during_Rie_subgiant.£hase_(Deheuvels 
et al. l20l4 hereafter DI4) and that a significant radial dif¬ 
ferential rotation establi shes by the time stars reach the 
red g iant branch (R GB) (iB eck et a ni2012l . iDeheuvels et al.l 
I20l4 hereafter DI2. lBeck et al. I2014D . This was expected 
as a consequence of the severe core contraction that occurs 
after the main-sequence turnoff. However, the core rota¬ 
tion rates measured for red giants are much slower than 
they would be in the absence of AM transport. Besides, 
by me asuring the core rota tion of several hundreds of red 
giants, iMosser et al.l (|2012al ) showed that even though stel¬ 
lar cores keep contracting, they actually spin down during 
the evolution along the RGB. These observations clearly 
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indicate that a very efficient mechanism of AM transport 
operates between the core and the envelope of red giants. 

For now it seems that none of the known mechanisms 
of AM transport are able to account on their own for the 
evolution of the rotation profiles of low-mass stars in the 
RGB obtained from seismology. Purely hydrodynamical 
processes including the effects of meridional circulation and 
turbulence shear as currently understood predict core rota¬ 
tion rates that ar e 2 to 3 orders of magnitude higher than 
the observations (lEggenberger et Eil]l2012l . iMaroiies et al.l 
1201,11 . ICeillier et all 201,lll . Magnetic instabilities were also 
considered as a possible source of AM transport. The insta¬ 
bility of toroidal magnetic fields under the influence of dif¬ 
ferential rotation (azimuthal mag netorotational instabili ty, 
AMRI) is a promising candidate. iRiidieer et al.l (l2014ai ra) 
showed that the resulting effective viscosity may be high 
enough to account for the slow core rotation rate of red 
giants, but stratification needs to be taken into account in 
order to.^uanti^^ts_effiect^_Oi^^he_othe£jiaiid^Cantiello 
et al. (I20l4i s howed that even if the so-called Tavler-Spruit 
(ISpruid I 2 OO 2 I ') dynamo loop exists, which is still debated, 
it does not generate sufficient AM transport to reproduce 
the slow core rotation rates observed in Kepler red giants. 
Interestingly, the transport of AM through IGW excited at 
the bottom of the convective envelope could account for the 
core/envelope decoupling during the subgiant phase. Dur¬ 
ing this phase, IGW are expected to damp out and release 
their AM before they reach the stellar core owing to the 
large increase in the Brunt-Vaisala fr equency, which could 
decouple the core from t he envelope (|Talon fc Charbonnell 
I 2 OO 8 L IfViller et aLll2014r ). However, IGW seem unable to 
then produce the recoupling that is necessa ry to account 
for the spin down of the core observed by iMosser et al.l 
(1201 2ah . 


So far, most of the constraints obtained from the seis¬ 
mology of red giants involve low-mass stars (M < 2 . 4 M 0 ) 
which ascend the RGB before triggering He-burning in a 
flash due to the degeneracy of the core. These stars are 
known as red clump stars in reference to their location in 
the HR diagram. Less constraints exist on th e internal ro¬ 
tation of intermediate-mass stars. Recently, iKiirtz et al.l 
(|2014r ) were able to seismically measure the internal rota¬ 
tion of a main sequence A-type star and found that this 
star rotates almost rigidly. However, it should be noted 
that this star might not be representative of all A-type 
stars since it has a peculiarly slow rotation rate. After the 
end of the MS, intermediate-mass stars cross the subgiant 
phase and the bottom of the RGB on a Kelvin-Helmoltz 
timescale, before settling in the clump (they are then re¬ 
ferred to as secondary clump stars). This raises the ques¬ 
tion whether AM transport can be fast enough to operate 
on such a short timescale. If not these stars should reach 
the clump with strong differential ro tation and very fast¬ 
spinning cores. iMosser et al.l (|2012af ) were able to seismi¬ 
cally measure the mean rotation rate in the g-mode cavity 
for secondary clump stars and found periods ranging from 
20 to 250 days. One must keep in mind that the He-burning 
core is convective, so gravity waves do not propagate in¬ 
side the_corean^^he,^otatio^^ates_obtaine^^^Mosser 
et al. (|20I2al l correspond to the layers just above the He- 
burning convective core. Assuming solid-body (SB) rota¬ 
tion a t the end of the main sequence. iTavar &: Pinsonneault 
( 201J|) showed that the rotation rates found bv IMosser et ah 
(20123) are compatible with a SB-rotation profile in the 


clump for intermediate-mass stars, which suggests that an 
efficient AM redistribution has operated by the time these 
stars reach the clump. 

In this paper, we aim at testing the hypothesis of a SB- 
rotation in the interior of intermediate-mass clump stars 
with seismology, using Kepler data. It has been shown that 
a ratio between the mean core rotation and the mean enve¬ 
lope rotation can be estimated from the splitt ings of mixed 
modes in subgiants and young red giants, see lGoupil et al.l 
(l201.ir i (hereafter G13), D14. On the other hand, for more 
evolved red giants the rotational splittings are dominated 
by the contribution of the central layers, and no reliable 
estimate of the envelope rotation can be obtained (G13). 
We here show that constraints on the amount of differen¬ 
tial rotation can be obtained for clump stars, provided the 
splittings of p-dominated modes can be reliably estimated. 
This latter condition is not obviously satisfied for clump 
stars because the mode widths of their p-dominated modes 
are comparable to the rotational splittings. We thus dedi¬ 
cate special care to establishing the significance of our esti¬ 
mates of the splittings of these modes. In Sect, dj we select 
the red-clump stars that are the most likely to provide con¬ 
straints on the amount of differential rotation among the 
Kepler targets. We then proceed to estimate the mode fre¬ 
quencies and rotational splittings in Sect. [2l and we address 
the question of the significance of the estimated splittings. 
In Sect, m we validate the simplified approach proposed 
by G13 to measure the internal rotation using the observed 
splittings and we apply it to all the stars of our sample. 
We thus find evidence for a weak radial differential rota¬ 
tion in six of the seven selected targets, the last one being 
consistent with a SB rotation. In Sect. Owe try to obtain 
complementary information on the rotation profiles of three 
stars of the sample based on high-resolution spectroscopic 
observations (HERMES spectrograph), but the constraints 
we obtain are too loose to test the seismic results. We fi¬ 
nally discuss the implications of our results in terms of AM 
transport in intermediate-mass stars in Sect. O 

2. Selection and characterization of targets 

Among the red giants observed with Kepler, we searched 
for secondary clump stars that are most likely to provide 
constraints on the amount of differential rotation. For this 
purpose, we used the complete Kepler dataset available, 
i.e. 4 years of data from quarters QO through Q17 with 
the long cadence mode (integration time of 29.4 min). The 
light cu rves were processed u sing the Kepler pipeline devel¬ 
oped bv iJenkins et al.l (|2010r) . and corrections from outliers, 
occas ional jumps, and drifts were also applied (|Garcfa et al.l 
120111 ). The power density spec tra were obtained by usiii g 
the Lomb-Scargle periodogram (|Lomblll97a . lScargldll982D . 
and we then applied the screening process that is described 
in the sections below in order to identify optimal targets. 


2.1. Selection based on stellar mass 

We first applied a mass criterion to select secondary clump 
stars. Standard stellar models predict that stars that form 
with a mass over 2.4 Mq trigger He-burning in a non¬ 
degenerate core. However, non-standard physical processes 
such as core overshooting during the M S can have a large 
influe nce on this theoretical threshold (|Montalban et all 
I 2 OI 3 D . Besides, stars are expected to lose mass during 
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their ascent of the RGB, so that the lower limit on the 
masses of secondary clump stars is in fact lower than the 
predicted theoretical threshold. Based on the properties 
of the observed mixed mod es in Kepler secondary clump 
giants, iMosser et al.l (|2014f ) observationally estimated the 
mass-limit to about 1.9 ± 0.2 Mq (see Fig. 1 of their pa¬ 
per). In this study, we used the upper bound of this obser¬ 
vational threshold in order to ensure that only secondary 
clump stars are retained, and selected stars with masses 
above 2.1 Mq. For this purpose, we applied seismic scal¬ 
ing relations, which required us to estimate the mean large 
separation of acoustic modes Av and the frequency of max¬ 
imum power of solar-like oscillations t'max in the considered 
targets. An estimate of the effective temperature is also 
needed, which we obtain photometrically following ? (see 
Sect. ITOl) . 
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2.1.1. Estimate of i^max 

In order to estimate the frequency of maximum power of 
the oscillations, we fitted a Gaussian envelope to the power 
spectrum in the neighborhood of the detected solar-like 
oscillations using a maximum likelihood estimate (MLE) 
method. The central frequency of the fitted Gaussian func¬ 
tion provides an estimate of t'max- The fit required to in¬ 
clude a description of the background. The contribution 
fro m granulation was modeled as a Harvev profile f Har- 
vev 119851 ). and white noise was a dded to account fm pho- 
ton noise. As wa.s poi nted out bv iKaroff et al.l (|20ll) and 
iKallinger et al.l ()2014l ). we found that an additional Har¬ 
vey prohle was required to correctly fit the observed back¬ 
ground. Beside providing an estimate of i^max, this fit yields 
a model for the background, which we further use to extract 
the properties of the oscillation modes (Sect. E]). 


2.1.2. Estimate of Aiy 


To estimate the mean larg e separation A i / of th e targets, we 
followed the procedure of iMosser et al.l ([20l3). It consists 
in htting a second-order asymptotic expression to the ob¬ 
served modes. In stead of i nclud ing the second-order terms 
as prescribed bv iTassoiiil (|l98r)l ). the authors proposed to 
introduce a constant parameter a corresponding to a linear 
variation of the large separation with radial order n. This 
yields the following expression for radial modes 


Vn,l=G = 


T “t” 77 ^max) 


Av 


( 1 ) 


The index Umax was chosen so that it matches the frequency 
of maximum power of the oscillations i^max) and it is there¬ 
fore not necessarily an integer. By smoothing the power 
spectrum of the observed stars, we obtained hrst approx¬ 
imate estimates of the frequencies of ^ = 0 modes. The 
expression given by Eq. [l]was then htted to our estimates 
of the I = 0 mode frequencies, yielding optimal values of 
Ah', £p, and a. An example of this fit is shown in Fig. 
[T] for one star of the sample. This provides a precise ob¬ 
servational value of Ah', which can be translated in the 
asymptotic value if needed. 


2.1.3. Seismic scaling relations 

We then used seismic scaling relations that relate the global 
seismic parameters Ah' and to stellar properties such 


Fig. 1. Echelle diagram of KIC7581399, obtained from 4 
years of Kepler data. The white diamonds correspond to a fit of 
the asymptotic expression of radial modes given by Eq. [T]to the 
observed radial modes. The gray and red diamonds correspond 
to a fit of the asymptotic expression of mixed modes given by 
Eq. [21 m and [S] (only the red modes were used to perform the 
fit). 


as the mass, radius and surface gravity (iBrown et al.lll99lT l. 
These relations are derived from the fact that i/max scales 
to good approxiination as the acoustic cut-off frequenc y 

(|Brown et al.lll99ll . IStello et 70120081 . iBelkacem et al.ll201ll ) . 

To obtain estimates of the stellar masses and radii of the 
considered stars, an estimate of the effective temperature 
of the star was also needed. We have used the photomet- 
ric estimates of Tpf f obtai ned from the recipe proposed by 
iPinsonneault et all ()2012D that we applied to the griz pho¬ 
tometry in the Kepler input catalog (KIC). As mentioned 
above, we retained only the stars that have a mass above 
2.1 Mq. 

2.2. Selection based on the evolutionary stage 

To select secondary clump stars, we needed to distinguish 
core-He-burning (clump) stars from H-shell-burning (RGB) 
stars. As is well known, this cannot be achieved merely 
from the location of stars in the HR diagram. However, 
these two classes of stars have a very different core struc¬ 
ture (clump stars have a convective core due to core-He 
burning). As a result, the mean period spacing of g modes 
for these two classes of stars are s ufficiently different t o 
distinguish them spectacularly well ([Bedding et al.l[20Tll ). 
We thus evaluated the period spacing AHi of Z = 1 modes 
in the cons idered stars. For thi s purpose, we followed the 
recipe that iMosser et al.l (|2012bl ) adapt ed from the asymp - 
totic expression of p-g mix ed modes of lUnno et ahl (|1989l ). 
According to Eq. 16.50 of lUnno et al.l (|l989l ). the match¬ 
ing in the evanescent zone of solutions corresponding to g 
modes in the core and to p modes in the envelope requires 
that 

tan(6(p) = gtan(0g) (2) 

where 

l-rb l-rd 

0g= kr dr, and 9p= kr dr (3) 


Article number, page 3 of llbl 





















































Table 1. Period spacings 


KIC-Id 


KIC8962923 

KIC5184199 

KIC4659821 

KIC3744681 

KIC7467630 

KIC9346602 

KIC7581399 


AIIi (s) q 


299.8 

9.27 

239.7 

9.22 

214.6 

9.18 

233.1 

9.41 

293.5 

9.19 

241.4 

9.23 

222.4 

9.18 



doi (fiH: 

9.76 

0.34 

9.91 

0.48 

9.66 

0.51 

9.45 

0.16 

9.85 

0.39 

9.43 

0.47 

9.74 

0.47 


and the radii ra and rb (resp. and rd) are the inner and 
outer turning points of the g-mode (resp. p-mode) cavity. 
The parameter q corresponds to the coupling between the 
p- and g- mode cavities. Its va riations with frequency are 
neglected. iMosser et al.l (|2012bf ) have developed the expres¬ 
sions of the phases 0-p and 0g in order to take into account 
higher-order terms in the asymptotic developments of p and 
g modes, yielding 


L 

0„ = TV 


^ - K)n.l 
1 


AIIii/ 


— £. 


(4) 

(5) 


where the ^ correspond to the frequencies of theo¬ 

retical pure Z = 1 p modes. The values of (t'p)„ ^ can be 
obtained from the second-order asymptotic expression of 
radial p modes as follows 


= ’^ri,i=o + (1/2 - (ioi)Ai/ (6) 

where doi corresponds to the mean small separation built 
with I = 0 and I = 1 modes. The value of this quantity 
is a priori not known since pure Z = 1 p modes cannot be 
observed. 

For a given period spacing AIIi, coupling term q, phase 
term Eg, and doi parameter, the mode frequencies can 
be approximated by the solutions of Eq. [5] that can be 
obtained with a Newton-Raphson algorithm. We then 
searched for the combination (AIIi, g, Eg, cZoi) that repro¬ 
duces the observed modes at best. The results are given 
for each star in Table [TJ The mode frequencies obtained 
with the optimal set of parameters for one target of the 
sample are shown in Fig. [TJ We note that we have fur¬ 
ther used our estimates of AIIi, q, and Eg to evaluate the 
trapping of the modes in Sect. 14.21 

We retained only the stars for which the mean period 
spacing AIIi was foun d to be above 120 s iii order to select 
only clump stars (see [Bedding et all 1201 ll . iMosser et all 
[ 2011 . 


2.3. Selection based on the rotational splittings of the modes 

Finally, an obvious condition to measure the amount of dif¬ 
ferential rotation was that the oscillation modes of the stars 
should be clearly split by rotation so that reliable rotational 
splittings can be extracted from the power spectra. A first 
inspection of Kepler oscillation spectra immediately shows 
that such stars are harder to find among clump stars than 
among RGB stars. There are two reasons for this: (1) the 
lifetimes of the modes are shorter for clump stars, so that 


they ha ve larger widths, and (2) as shown bv iMosser et al.l 
(l2012bf l. the cores of clump stars have longer rotation peri¬ 
ods than those of RGB stars, which means that their rota¬ 
tional splittings are smaller. As a result of these two com¬ 
bined effects, the rotational splittings of the modes are often 
comparable to their widths, which complicates the mea¬ 
surement of splittings. This is not so much a problem for 
g-dominated modes, which have larger inertias and there¬ 
fore shorter widths. However, the situation is more compli¬ 
cated for p-dominated modes. We stress that the splittings 
of these latter modes are absolutely necessary to estimate 
the amount of differential rotation between the core and 
the envelope. We therefore performed a first pre-selection 
of targets based on a visual inspection of the echelle dia¬ 
grams of all the considered giants to identify those whose 
mixed modes are most clearly split. We note that the most 
favorable stars are those that are seen equator-on (inclina¬ 
tion angle i = 90°), because in this case only the to = ±1 
components of Z = 1 modes, which are separated by twice 
the rotational splitting, are visible. 

The characteristics of the selected targets are listed in 
Table HI Only the targets that provided enough reliable 
splittings (see Sect. |3|) to estimate differential rotation are 
mentioned. The following study is based on this sample of 
seven targets. 


3. Extraction of rotational splittings 

In order to measure the amount of differential rotation in 
the selected targets, we needed to extract the rotational 
splittings of the modes. For this purpose, we used a max¬ 
imum likelihood estimate (MLE) method in the same way 
as D14. We here only briefly summarize it. We performed 
individual fits of each rotational multiplet (modes of same 
radial order n and degree Z). Among the multiplets, the to- 
components were modeled as Lorentzian profiles. The com¬ 
ponents were assumed to have a common width F. Within 
the multiplets, the height ratios between the components 
are given by a visibility factor that depends on the inclina¬ 
tion a ngle i of the star (|Gizon fc Solankill2003l . [Ballot et al.l 
1299611 . Finally, we assumed that the multiplets are sym¬ 
metric with respect to the central to = 9 component. They 
are thus equally spaced by the rotational splitting whose 
expression is given by 

= (7) 

where the functions Kn,i(r) are the rotational kernels, 
which depend on the eigenfunctions of the modes. 

As was mentioned in Sect. HJ one of the main chal¬ 
lenges of this work was to establish that the splittings of 
p-dominated modes can be recovered for clump stars. For 
the selected targets, several p-dominated modes are visually 
split (see Fig. HjQ. However, owing to the larger linewidths 
of these modes, it is not clear at first sight whether this is 
caused by the rotational splitting of the modes, or only the 
result of stochastic excitation. To distinguish between these 
two cases, we performed statistical tests, following both a 
frequentist and a Bayesian approach in a complementary 
way. 

^ The trapping of the modes is evaluated in Sect. |T] 
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Table 2. Global parameters of selected targets 


KIC-Id Al^ (^Hz) Teff (K) M/Mq R/Rq logs 


KIC8962923 

KIC5184199 

KIC4659821 

KIC3744681 

KIC7467630 

KIC9346602 

KIC7581399 


6.95 ±0.13 
7.85 ±0.13 
8.19 ±0.09 
5.47 ±0.13 
6.06 ±0.09 
5.07±0.13 
6.57 ±0.13 


79.3 ±0.3 
94.7 ±0.4 
101.5 ±0.3 
61.1 ±0.5 

70.6 ±0.4 

55.7 ±0.3 
80.9 ±0.4 


5253 ±91 
5197 ±89 
5160 ±98 
5084 ± 68 
5114 ± 103 
5162 ± 79 
5301 ± 80 


2.17 ±0.25 
2.23 ±0.23 
2.29 ±0.19 
2.46 ±0.35 
2.55 ±0.27 
2.58 ±0.37 
2.92 ±0.35 


9.36 ±0.48 
8.70 ±0.40 
8.53 ±0.30 
11.45 ±0.73 
10.81 ±0.49 
12.22 ±0.80 
10.73 ±0.57 


2.832 ±0.013 
2.907 ±0.012 
2.935 ±0.013 
2.712 ±0.015 
2.776 ±0.015 
2.675 ±0.013 
2.843 ±0.013 


3.1. Frequentist approach 

The frequentist approach consists in computing a false- 
alarm probability using the so-called Hq test. For each 
of the observed I = 1 multiplets in the oscillation spectra 
of the selected targets, we performed two MLE fits using 
successively the two following hypotheses: 

— Hq hypothesis: the observed multiplet has a splitting 
much smaller than the mode width {Sv^ <C F). In 
this case, the mean profile of the whole multiplet cor¬ 
responds to a single Lorentzian prohle and is described 
by no = 3 parameters (frequency vq of the central com¬ 
ponent, linewidth F, and height H). 

— Hi hypothesis: the rotational splitting of the mode is 
at least of the order of magnitude of the mode width 
{Si's > F). The mean profile then corresponds to a well 
separated multiplet and is described by ni = 5 parame¬ 
ters (same ones as for the Hq hypothesis plus the incli¬ 
nation angle i and the splitting Svs). 

Let us denote as io and £i the likelihoods of the best fits 
obtained under the Hq to the Hi hypotheses, respectively. 
We note that fitting the power spectrum in the neighbor¬ 
hood of the considered multiplet under the Hi hypothesis 
necessarily gives an agreement with the data at least as 
good as fitting under the Hq hypothesis. Indeed, the best 
ht under the Hq hypothesis can also be obtained with the 
Hi hypothesis and a pole-on inclination (i = 0°). Asa con¬ 
sequence, the likelihood £o is always lower than the likeli¬ 
hood ii. To determine whether or not the Hq hypothesis 
can be rejected, one needs to assess the significance of the 
likelihood impr ovemen t when switching from the Hq to the 
Hi hypothesis. IWilksI (Il938f ) has shown that the quantity 
AA = —2 [ln(.^l) — ln(£o)] follows the distribution of a 
with (ni — Tin) degrees of freedom fsee also Appourchaux 
et al. 119981 1. We have verified in our particular case that 
AA indeed follows the distribution of a with ni—riQ = 2 
degrees of freedom by performing Monte Carlo simulations 
(see Appendix 121) . 

Supposing that we have computed the value AA for 
a given mode, the false-alarm probability is estimated by 
computing the probability that a mode in the Hq hypoth¬ 
esis (i.e. not rotationally split) can produce such a high 
value of AA. This probability is given by 

P = -P(x^(2 degrees of freedom) > AA) (8) 

and is commonly referred to as the p-value. Low p-values 
indicate that a mode described as a sing Lorentzian profile 
is unlikely to produce the observations. 


For each observed I = 1 mode, we computed the p-value 
using Eq. [51 and we rejected all the modes for which a 
p-value higher than 0.05 was found. We applied this pro¬ 
cess to all the targets that were pre-selected in Sect. [H 
and kept only the stars for which at least hve signihcant 
splittings were obtained. This left us with seven targets, 
which are listed in Table [H with a number of signihcant 
splittings ranging from 6 to 14. The extracted frequen¬ 
cies and rotational splittings of the selected multiplets are 
given along with their p-values in Tables [31 and [H for all 
the stars of this sample. To each multiplet corresponds an 
estimate of the inclination angle i. We note that for ev¬ 
ery star the angles found for each mode agree quite well 
with one another and consistently point toward inclination 
angles close to 90° (see Tables [3] and SI). This is under¬ 
standable considering that the stars were selected so that 
the modes are clearly split, which greatly favors stars with 
high inclination angles, as pointed out in Sect. S) We note 
that in two cases (mode at v ^ 89.1 p,Hz for KIC7581399 
and mode at ^ ~ 92.6p,Hz for KIC5184199), the best-ht 
solutions point toward a solution with an inclination an¬ 
gle around i = 55° and a splitting around twice the value 
obtained for the other modes. This is likely caused by the 
influence of a higher-degree mixed mode nearby (see discus¬ 
sion about these modes in Sect. [33), and these two modes 
were discarded in the following. 

3.2. Bayesian approach 

The limits of the Hq test used in the frequentist approach 
are well known. To summarize in our case, a low p-value en¬ 
sures that noise (in our case a multiplet where the splitting 
is much smaller than the mode width) has little chances 
to produce the observed signal. However, it says nothing 
about the probability to produce the observations under 
the Hi hypothesis (i.e. when the multiplet is detectably 
split by rotation). 

This can be remedied by applying a Bayesian approach, 
which produces an estimate of the consistency between a 
given model Mj and a dataset D. In the present case, the 
dataset D corresponds to the observed power spectrum in 
the neighborhood of a given detected I = 1 mode, and we 
aim at testing the two models of power spectrum Mq and 
Ml that correspond to the hypotheses Hq and Hi defined 
in Sect o In other words, model Mq is built assuming of 
a single Lorentzian prohle, whereas model Mi corresponds 
to a triplet of Lorentzian prohles. Each model Mj relies on 
a set of free parameters 0, which contains for both models 
the background level, the central frequency t'o, the mode 
linewidth F, and the mode amplitude A, and additionally 
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Fig. 2. Section of the oscillation spectrum of KIC3744681 in the neighborhood of 4 p-dominated modes. The red curves 
represent the best fits to the data. 


for model Mi the inclination angle i and the rotational 
splitting 8vs. 

Within a Bayesian approach, we must define a prior 
probability law TT[6\Mj, I), based on a certain number of 
assumptions on the signal that are here summarized as I. 
The chosen priors are explicited below. We used the Bayes 
theorem to define the posterior probability density function 
of a given set of parameters 6 assuming model Mj and 
knowing the dataset D through the relation 


7r(0|Z?,M„J) 




(9) 


where 'n{D\6,Mj,I) corresponds to the likelihood already 
introduced in Sect. [SIl ■n{9\Mj,I) corresponds to the 
afore-mentioned priors, and P{D\M, I) is the evidence of 
model Mj knowing the data D and is defined as 

P{D\Mj,I)= J T:{e\Mj,I)T:{D\e,Mj,I)de. (10) 

We assumed the following prior probability laws for the 
free parameters 9: 

— amplitude A: uniform prior for the variable over 

[ln(Po/1000), ln(1.2Po)] where Pq is the total power in 
the fitting window. This non-informative prior is very 
loose. 

— width T: uniform prior for InT over the range 
[Inr, ln(0.4 p.Hz)] with Gaussian decays around these 
boundaries, r is the frequency resolution of the spec¬ 
trum (about 0.008 /iHz) and 0.4 /rHz corresponds to 
an upper limit of radial mode widths observed in these 
stars. With this prior, we assumed that the widths of 
I = 1 mixed modes cannot excess those of Z = 0 modes, 
which are purely acoustic. 


— central frequency vq: uniform prior for vq spanning an 
interval of 0.6 /rHz around the middle of the fitting win¬ 
dow. 

— rotation splitting Sv^: uniform prior for Svs over 
[0,0.3 ^Hz]. 

— inclination angle i: uniform prior for cosi over [0,1]. 
This prior derives from the assumption of isotropy of 
the inclinations of stars. 

— background level: we used as a prior the result of a 
global fit of the background as described in Sect. 13.11 
By doing so, we used information coming from outside 
the fitting window. 

To compare two different models Mi and Mj we used 
the evidence of these models as introduced in Eq. (TU] to 
define the so-called odd ratio 

P{M,\D,r) P{M,\I) P{D\M,,I) 

P{M,\D,I) PiM,\I)P{D\M„iy ^ 

By assuming that the two competing models Mi and Mj 
are equiprobable, i.e. P{Mi\I) = P{Mj\I) = 1/2, the odd 
ratio is equivalent to the Bayes factor 


, _ p{Dm,i) 

P{D\Mj,I)' 


( 12 ) 


The relative probability to favor model Mj over Mi is then 
given as 

P{MjD) = {l + 0.y-^ (13) 

which derives from Eq. [TT]if we assume that no other model 
than models Mi and Mj is possible. 
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Table 3. Extracted frequencies uq and rotational splittings Sus for the 1 — 1 modes of K1C8962923, KIC5184199, KIC4659821, 
and KIC3744681. The 3'^'^ column gives the estimate of the inclination angle i obtained from each mode. Only modes whose 
splitting was found significant with a p-value (given in 4*** column) below 0.05 are given. The S*** column gives the posterior 
probability obtained from a Bayesian analysis that the rotational splitting of the mode is detected (see Sect. EH). The modes 
whose parameters are shown in italics were discarded (see text). 


KIC8962923 


r'o (mHz) 

Si's (mHz) 

*(°) 

p-value 

Posterior probability (%) 

57.275 ±0.004 

0.070 ±0.004 

70.68 ±8.97 

0.015 

49.3 

70.532 ±0.004 

0.069 ±0.004 

81.46 ±6.65 

0.000 

99.1 

71.883 ±0.006 

0.058 ±0.006 

90.00 ± 14.54 

0.004 

98.3 

72.932 ±0.010 

0.080 ±0.010 

90.00 ± 10.66 

0.027 

91.2 

73.958 ±0.004 

0.065 ±0.005 

68.92 ±6.47 

0.001 

99.5 

77.046 ±0.004 

0.077 ±0.004 

72.66 ±6.30 

0.005 

97.5 

78.672 ±0.006 

0.083 ±0.006 

65.72 ±7.83 

0.020 

94.9 

81.110 ±0.008 

0.056 ±0.008 

90.00 ± 15.31 

0.042 

67.9 

84.792 ±0.007 

0.079 ±0.007 

80.62 ±7.29 

0.001 

99.8 


KIC5184199 


uo (mHz) 

Si's (mHz) 

* (°) 

p-value 

Posterior probability (%) 

69.466 ±0.009 

0.105 ±0.009 

90.00 ±9.60 

0.019 

92.0 

76.319 ±0.008 

0.086 ±0.009 

90.00 ±6.70 

0.003 

95.0 

82.640 ±0.011 

0.098 ±0.011 

90.00 ± 12.98 

0.039 

55.8 

84.883 ±0.009 

0.080 ±0.009 

89.99 ±67.88 

0.008 

99.1 

91.505 ±0.007 

0.093 ±0.007 

90.00 ± 10.56 

0.000 

100.0 

92.614 ± 0.009 

0.130 ± 0.012 

54.09 ± 6.90 

0.010 

86.2 

94.115 ±0.010 

0.102 ±0.009 

90.00 ±9.93 

0.000 

98.2 

98.241 ±0.004 

0.093 ±0.004 

79.98 ±4.10 

0.000 

100.0 

101.195 ±0.009 

0.076 ±0.010 

90.00 ± 17.27 

0.003 

96.7 

105.753 ±0.014 

0.106 ±0.013 

90.00 ± 17.29 

0.043 

99.6 


KIC4659821 


r'o (mHz) 

Svs (mHz) 

i(°) 

p-value 

Posterior probability (%) 

80.659 ±0.007 

0.088 ±0.008 

74.27 ±7.66 

0.002 

93.2 

94.645 ± 0.004 

0.082 ±0.005 

78.71 ±5.16 

0.000 

99.8 

96.193 ±0.009 

0.069 ±0.009 

90.00 ± 16.60 

0.019 

91.7 

110.366 ±0.008 

0.085 ±0.008 

90.00 ±9.88 

0.006 

99.5 

112.513± 0.012 

0.082 ±0.017 

79.58 ±22.04 

0.011 

83.5 

113.765 ±0.010 

0.084 ±0.013 

71.38 ± 10.90 

0.006 

94.8 


KIC3744681 


Vo (mHz) 

Svs (mHz) 

i n 

p-value 

Posterior probability (%) 

51.747 ±0.002 

0.093 ±0.002 

70.30 ±8.08 

0.000 

71.6 

52.286 ±0.009 

0.108 ±0.010 

77.16 ±6.49 

0.018 

75.5 

57.558 ±0.010 

0.082 ±0.009 

90.00 ±8.51 

0.015 

92.6 

58.093 ±0.007 

0.082 ±0.007 

90.00 ± 17.37 

0.000 

96.4 

61.856 ±0.004 

0.090 ±0.004 

73.38 ±9.17 

0.000 

99.3 

62.822 ±0.007 

0.095 ±0.007 

79.48 ±5.76 

0.000 

99.7 

63.542 ±0.009 

0.088 ±0.010 

74.02 ±8.09 

0.000 

97.8 

64.215 ±0.009 

0.109 ±0.009 

80.26 ±8.43 

0.001 

93.3 

67.845 ±0.009 

0.097 ±0.009 

90.00 ± 142.99 

0.000 

98.8 

68.803 ±0.010 

0.074 ±0.011 

80.16 ± 12.18 

0.017 

70.8 

69.560 ±0.011 

0.121 ±0.011 

69.23 ±7.10 

0.020 

72.6 

73.660 ±0.008 

0.088 ±0.007 

90.00 ± 11.93 

0.000 

92.3 

74.545 ±0.009 

0.100 ±0.010 

67.17 ±7.29 

0.004 

58.8 

75.537 ±0.013 

0.098 ±0.011 

79.17 ± 12.22 

0.036 

69.1 
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Table 4. Same as Table |3] for targets KIC7467630, KIC9346602, and KIC7581399. 


KIC7467630 


Vo (pHz) 

Svs (pHz) 


p-value 

Posterior probability (%) 

61.557 ±0.002 

0.059 ±0.001 

58.65 ±47.98 

0.032 

60.7 

67.666 ±0.004 

0.064 ±0.003 

89.99 ± 16.34 

0.000 

59.0 

68.903 ±0.008 

0.064 ±0.007 

90.00 ± 13.72 

0.011 

97.1 

69.891 ±0.010 

0.076 ±0.010 

89.99 ±23.16 

0.005 

94.6 

70.769 ±0.008 

0.071 ±0.008 

90.00 ±20.24 

0.003 

98.7 

82.877 ±0.010 

0.101 ±0.012 

69.44 ±7.59 

0.011 

70.3 


KIC9346602 


Vo (pHz) 

SVa (pHz) 

in 

p-value 

Posterior probability (%) 

43.638 ±0.005 

0.073 ±0.005 

82.85 ±7.72 

0.000 

96.4 

44.324 ±0.005 

0.067 ±0.005 

90.00 ± 10.53 

0.000 

98.6 

48.203 ±0.006 

0.083 ±0.006 

90.00 ± 14.14 

0.001 

91.8 

48.688 ±0.007 

0.070 ±0.007 

78.18 ±8.28 

0.007 

70.2 

53.694 ±0.004 

0.080 ±0.004 

90.00 ±4.77 

0.000 

100.0 

54.119 ±0.005 

0.074 ±0.005 

82.59 ±6.82 

0.000 

99.8 

54.678 ±0.004 

0.077 ±0.005 

78.11 ±5.24 

0.000 

99.4 

57.587 ±0.007 

0.098 ±0.008 

76.24 ±6.34 

0.023 

78.7 

58.334 ±0.007 

0.073 ±0.007 

90.00 ±5.94 

0.006 

91.6 

60.037 ±0.003 

0.071 ±0.004 

77.90 ±5.35 

0.019 

58.1 

62.775 ±0.003 

0.080 ±0.003 

72.70 ±7.46 

0.000 

99.4 

64.942 ± 0.008 

0.074 ±0.008 

90.00 ± 17.27 

0.038 

65.6 

67.956 ±0.005 

0.082 ±0.005 

81.02 ±7.68 

0.000 

90.1 

73.922 ±0.002 

0.079 ±0.003 

78.09 ±6.60 

0.000 

97.5 


KIC7581399 


Vo (pHz) 

Svg (pHz) 

in 

p-value 

Posterior probability (%) 

63.546 ±0.006 

0.084 ±0.006 

74.69 ±7.28 

0.005 

95.7 

67.852 ±0.006 

0.086 ±0.004 

90.00 ± 10.37 

0.001 

50.5 

68.811 ±0.005 

0.079 ±0.005 

80.13 ±8.35 

0.031 

49.7 

76.375 ±0.010 

0.083 ±0.011 

78.07 ± 12.43 

0.003 

88.8 

81.297 ±0.005 

0.082 ±0.005 

90.00 ± 16.03 

0.000 

92.3 

82.502 ±0.009 

0.072 ±0.012 

78.81 ± 15.49 

0.023 

89.4 

83.314 ±0.010 

0.095 ±0.009 

89.98 ± 19.56 

0.001 

98.4 

89.056 ± 0.010 

0.155 ± 0.011 

60.01 ± 6.93 

0.001 

94.7 


To calculate the evidences of models Mq and Mi, we 
developed an automated parallel tempering Markov-Chain 
Monte Carlo similar to the one developed bv iBenomar et all 
(120091 1. We then computed the odd ratio Oqi for all the 
modes that were found to have a p-value below 0.05 in Sect, 
o and derived the probability P{Mi\D) to favor model 
Ml. The results are listed in the last column of Tables [3] 
and m All the modes for which the obtained probability 
was higher than 85% were considered as reliable and used 
to probe the rotation profiles of the selected targets. 

We stress that the approach followed here is conserva¬ 
tive. In particular, we have assumed that either model Mg 
or model Mi is correct, i.e. that the contribution from only 
one I = 1 mode is present in the htting window. If other 
non-negligible components (such as higher degree modes) 
are also present then model Mi can appear nearly as bad 
as model Mg in describing the spectrum, even though a 
clear I = 1 multiplet is present. This occurs especially for 


g-dominated modes, which lie in a region of the spectrum 
where mixed I = 2 and I = 3 modes can be present. For 
these modes, the relatively low posterior probability does 
not mean that the splitting cannot be measured, but that 
the model would need to be improved. In the present case, 
we did not seek to improve it because enough g-dominated 
modes were found signihcant with the current approach to 
estimate the amount of differential rotation (see Sect. U]). 

By imposing that the measured splittings should pro¬ 
duce both a low p-value (frequentist approach) and a high 
posterior probability P{Mi\D) (Bayesian approach), we en¬ 
sured that only very robust splitting detections are vali¬ 
dated, at the expense of potentially discarding good data. 
However, a robust selection process was needed to clearly 
establish that the splittings of both p and g-dominated 
modes can be reliably extracted, which is a necessary condi¬ 
tion to estimate the amount of differential rotation. With 
this double selection process, we retained between 3 and 
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10 significant splittings for the stars of the sample. In the 
following only these splittings are used to infer the inter¬ 
nal rotation profile, but the splittings that passed only the 
frequentist test are also shown on the plots. 

4. Measuring the amount of differential rotation 

Two different approaches have been proposed to obtain in¬ 
formation on the internal rotation profile of red giants using 
the mode splittings that were selected in Sect. |31 

1. The first one requires to obtain a stellar model that 
closely matches the observed mode frequencies, as well 
as the classical constraints on the star (atmospheric pa¬ 
rameters). If such a model can be found, one gains 
access to the eigenfunctions of the modes, which makes 
it possible to determine where the modes are trapped 
and to compute the rotational kernels of the modes. 
Information on the rotation profile can then be ob¬ 
tained either by performing inversions of Eg. [7] (D12. 
D14) , or by applying forward techniques (|Beck et al.l 
I 2 OI 4 I I. In spite of these advantages, this method is time- 
consuming because it requires to compute a model that 
matches the observed mode frequencies. More impor¬ 
tantly, it is model-dependent. This can be a problem 
here since the high values of AIIi that are observed for 
a good proportion of Kepler clump stars fail to be re- 
orod uced bv several evolution codes te.g. Montalban 
et al. l201,‘l^. although some co des might not have this 
problem ( Lagarde et al.l[2M^ . It was suggested that 
this discrepancy could be solved by adding overshoot¬ 
ing at the b oundary of the convecti ve core induced by 
He-burning (|Montalban et al']l2013r ). but this is still a 
matter of debate. 

2. A simplified approach was pro posed by G13. By using 
asymptotic developments from IUnno et al.l (|1989f ). they 
showed that estimates of the mean rotation rate in the 
g-mode and p-mode cavities can be obtained using only 
the mode frequencies, i.e. without having to compute 
a stellar model of the observed target. This constitutes 
an appealing alternative to the first method in the cases 
where the observed AIIi cannot be reproduced by stel¬ 
lar models. The downside is that it relies on several 
layers of approximations and needs to be tested in the 
case of clump stars. 

We applied both approaches to one test-target of our 
sample. We chose KIC7581399 because it has a value of 
AIIi that is low enough to be satisfactorily reproduced by 
current stellar models, so that the first approach can be 
applied. We note that the internal structure of the best-fit 
models obtained for KIC7581399 fsee ld.lll is probably not 
entirely correct, because we know that current models un¬ 
derestimate the values of AIIi of clump stars. However, it 
has been shown that the rotation profiles obtained from in¬ 
versions do not critically depend on the input physics of the 
reference model, provided the frequencies of the observed 
modes are well reproduced (D12, D14). Our objectives were 
twofold. First, we wanted to verify that the rotational split¬ 
tings that were extracted in Sect. [3] are sufficient to provide 
a constraint on the ratio between the core and the envelope 
rotation rates. Secondly, we aimed at testing the simplified 
method of G13 on red clump stars, in order to determine 
whether it can be applied to all targets. 



Fig. 3. Echelle diagram of the best-fit model for KIC7581399. 
The open circles correspond to the observed frequencies and the 
colored filled symbols represent those of the model (blue squares: 
I = 0 modes, red triangles: I = 1 modes). 


4.1. Model-dependent method: the test esse of KIC7581399 

4.1.1. Modeling KIC7581399 

In order to perform inversions of Eq. [71 one needs to have 
access to the rotational kernels of the modes, and there¬ 
fore to a stellar model of the sta r For this purpos e, we 
used the evolutionary code MESA (iPaxton et al.ll2011 ). We 
used the OPAL 2005 equation of state ([Rogers fc NavfonovI 
l2002f ) and opacity tables. The nuclear react ion rates were 
comp uted using the NACRE compilation (jAngiilo et al.l 
Il999h . We assumed the class i cal so lar mixture of heavy 
elements of iGrevesse fc No^ (|1993D . Convection was de¬ 
scribed using^he_^lassica^^ixing_Jengt^^heor^_(B6hm- 
Vitense Tl958f l with a mixing length parameter calibrated on 


the Sun (omlt = 2.08). The effects of microscopic diffusion 
and overshooting from convective regions were neglected in 
this study. 

To mo del KIC7581399, we f ollowe d the procedure pro¬ 
posed by iDeheuvels fc Michell (|201iri to model subgiants 
with mixed modes, and later adapted to red giants by D12. 
We refer to these papers for more details about the method, 
which is only briefly recalled here. The authors showed that 
the combined knowledge of the mean large separation (Az/) 
and mean period spacing of ^ = 1 modes (AHi) yields tight 
constraints on the stellar mass and age when other physical 
inputs are fixed. To model a red giant, they thus advocated 
to vary all physical parameters except for the mass and age, 
which are determined each time through an iterative pro¬ 
cess so that the observed {Ai/) and (AHi) are correctly 
reproduced. In practice, the constraint given by (Au) is 
advantageously replaced by the observed frequencies of the 
lowest-order radial modes, because the latter are less sen¬ 
sitive to near-surface effects. For the models, the mean 
period spacing (AHi) was estimated from the approximate 
expression 




(14) 


We note that (AHi) could in principle be estimated from 
the mode frequencies themselves by following the procedure 
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Fig. 4. Extracted rotational splittings for I — 1 modes of 
KIC7581399 (open black circles indicate modes that passed both 
frequentist and Bayesian significance tests and open gray circles 
modes that passed the frequentist test only). The colored cir¬ 
cles indicate theoretical splittings for linear rotation profiles that 
reproduce the rotation of KIC7581399 in the g-mode cavity as 
obtained from the OLA method, and with ratios flcore/fisurf of 
1 (purple), 2 (blue), 10 (red), or 100 (green) from top to bottom. 


described in Sect. \T7[ This would however require to com¬ 
pute the frequencies of Z = 1 modes for a large number of 
models, which is time-consuming. It was shown that there 
is a good agreement with the estimate of (AIIi) obtained 
from Eg. 1141 and tha t obtained from the mode frequencies 
(|Mosser et al.ll201^ . D14). 

We computed a grid of models, varying both the initial 
helium content 1) and the initial metallicity {Z/X)i. For 
each point of this grid, the procedure described above was 
applied to estimate the mass and age. The mode frequen¬ 
cies of each model were then com puted using the oscillation 
code Lose (|Scufiaire et abll^OOSl ). Before comparing them 
to the observed frequencies, we corrected them from the 
we ll-known near-s u rface effects using the recipe proposed 
bv iKieldsen et al.l (120081 ). which consists of adding to the 
mode frequencies a power law whos e exponent is ca li brate d 
in the Sun. As was pointed out bv iKieldsen et al.l ([20^, 
mixed modes have a larger inertia in the core than pure 
acoustic modes. Therefore, for these modes the correction 
terms were weighted by the ratio between the inertia of the 
closestj;adiaI^Qd§AQAhgJngrti^^f_|hejnodfijtse]f_{Aerts 
et al. l20inl) . 

We found a model that reproduces the surface param¬ 
eters of the stars and closely matches the observed mode 
frequencies and the star’s effective temperature. An echelle 
diagram of its mode frequencies corrected from near-surface 
effects as described above is shown in Fig. [31 This model, 
whose main properties are listed in Table [3] is used as a 
reference model in the following. 


4.1.2. Mode trapping 

The eigenfunctions of our reference model of KIC7581399 
were used to estimate the trapping of the modes. For each 
mode, we computed the parameter ( defined as the ratio 
between the kinetic energy of the mode in the g-mode cavity 


Fig. 5. Averaging kernel obtained for KIC7581399 with the 
OLA method in the g-mode cavity. 


and the total kinetic energy 




dr 


l(l 


1 )^] dr' 


(15) 


where and are the radial and horizontal displacements. 
A value of C close to 1 indicates that the mode is mainly 
trapped in the g-mode cavity, and a value of C close to 0, 
that it is trapped in the p-mode cavitj0. Fig. [^represents 
the rotational splittings of KIC7581399 that passed the sta¬ 
tistical test of significance described in Sect. [3|as a function 
of the mode trapping C. We clearly observe that the split¬ 
tings of p-dominated modes are comparable to the splittings 
of the g-dominated modes. This is a striking difference with 
the RGB stars, for which p-dominated mode s have a split¬ 
ting a bout half that of g-dominated modes (jMosser et al.l 
l2012al ). This already suggests that there might be less dif¬ 
ferential rotation in secondary clump stars than in RBG 
stars. 


4.1.3. Rotation in the g-mode cavity 

It was shown for RGB stars that the mean rotation rate 
in the g-mode cavity can be prec isely determined usin g the 
splittings of mixed modes (D12. [Mosser et n]l2012^ . For 
this purpose it is very convenient to use the OLA (optimally 
localized averages) method to invert Eq. [71 This method. 


et al. 19981 . iGhaolin et al.lll999l ) as well as to subgiants and 
RGB stars (D12, D14), consists of building combinations 
of the rotational kernels such that the resulting averaging 
kernels /C(r; ro) = Ck{ro)Kk{r) are as localized as possi¬ 
ble around a target point rg. When setting rp = Tg, where 
rg is the mean radius of the g-mode cavity, we obtain an 
averaging kernel that is well localized in the g-mode cavity 
(see Fig. 0. In particular, the contribution from the enve¬ 
lope is efficiently suppressed. This can be used to obtain a 
precise estimate of the mean rotation in the g-mode cavity 


^ The p-mode cavity depends on the degree I of the mode. Since 
we here consider 1 = 1 modes only, the p-mode cavity refers to 
the cavity of / = 1 p modes throughout the paper. 
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Table 5. Parameters of the best-fit model of K1C7581399 


Mass (Mq) 

Age (Myr) 

hi {Z/X)i 

Radius (Rq) 

TeS (K) 

2.65 

5686 

0.29 0.035 

10.38 

5129 


(rig). For KIC7581399, we obtained (ng)/(27r) = 173 ± 20 
nHz. 

We must stress that contrary to RGB stars, the core 
of clump stars is convective, due to He-burning. The con¬ 
vective core is very small because of the high temperature 
sensitivity of the 3a reaction (it represents less than 0.2% 
of the stellar radius and less than 4% of the stellar mass for 
KIC7581399) but nonetheless, the most central layers are 
not probed by g modes, which are evanescent in convective 
regions. So we must keep in mind that the rotation in the 
very core cannot be probed by seismology. 


4.1.4. Core-envelope contrast 

One crucial point in this study, which aims at estimating 
the amount of radial differential rotation in clump stars, 
is to determine whether or not we can also obtain an esti¬ 
mate of the mean envelope rotation rate. For this purpose, 
we considered simple "two-zone" models that rotate as solid 
bodies with a rate fig in the g-mode cavity, and with a rate 
rip in the p-mode cavity. Naturally, we used for fig the 
estimate that was obtained from the observations in Sect. 
14.1.31 We chose values of flp such that the ratio flg/flp 
is equal to 1, 2, 10, and 100 successively. These synthetic 
rotation profiles along with the rotational kernels of the 
best-fit model for K1C7581399 were plugged into Eq. [7] to 
obtain theoretical splittings, which are plotted as a func¬ 
tion of the mode trapping ( in Fig. U) This figure clearly 
shows that the splittings of p-dominated modes can be 
used to constrain the mean rotation in the p-mode cavities. 
For K1C7581399, a SB rotation profile throughout the star 
(Clg = Clp) can be excluded, as well as rotation profiles with 
a large amount differential rotation (flg/rip > 10). The 
observations appear to be consistent with a core-envelope 
rotation contrast around 2 (blue curve). 

This can be further quantified by searching for the pa¬ 
rameters (fig, rip) that reproduce the observed splittings 
at best. For this purpose, we minimized the reduced x^ed 
defined as 

1 ^ 

_2 X! 

fc=l 

where 6 vk are the M observed splittings, Cfe the corre¬ 
sponding error bars, rg is the limit between the p-mode 
and g-mode cavities, Ak{ro) = /J^“iFfc(r)dr and Bk{ro) = 

Kk{r) dr. We note that rg can be chosen anywhere in 
the evanescent zone between the p-mode and g-mode cavi¬ 
ties with very little change in the results. For KIC7581399, 
we obtained flg/( 27 r) = 156±12nHz and flp/( 27 r) = 95±15 
nHz. The mean rotation rate in the g-mode cavity is 
fully consistent with the result obtained with the OLA 
method. We thus obtained a core-envelope rotation con¬ 
trast of 1.6 ± 0.4 for KIC7581399. This shows that the 
amount of differential rotation in the secondary clump star 
K1C7581399 is much lower than for RGB stars. 


5vk - flpA/c(rg) - Vt^Bkiro) 




(16) 


Table 6. Mean rotation rates in the g- and p-mode cavities 
obtained for the selected targets with the method of G13 


KIG-Id 

(nH.) 


{%) 

8962923 

138 ±8 

79± 10 

1.8 ±0.3 

5184199 

200 ± 13 

63 ±20 

3.2 ± 1.0 

4659821 

165 ± 14 

79± 15 

2.1 ±0.4 

3744681 

194 ± 20 

63 ±36 

3.1 ±1.8 

7467630 

121 ±18 

96 ±28 

1.3 ±0.4 

9346602 

164 ±6 

53± 15 

3.1 ±0.9 

7581399 

167 ± 11 

84± 14 

2.0 ±0.4 


4.2. Testing the lGouoil et al\ K20li) approach on 
KIC7581399 

We then tested the simplified model-independent approach 
proposed by G13 on our test-target K1G7581399 in order 
to determine whether it can be applied to the targets of 
our sample. The method of G13 relies on two consecutive 
approximations. 

The first step consists in approximating the mode trap¬ 
ping parameter C by an expression ( that depends solely 
on the mode frequencies, with the advantage that C can be 
estimated from the observations. To establish this relation, 
the authors used approximate expressions of the eigenfunc¬ 
tions derived from JWKB analysis. For all the modes whose 
splittings could be reliably extracted in KIG7581399, we 
computed C following the prescription of G13 (Eq. A.27 
and A.28 of their paper). As can be seen in Eig. [S] (gray 
symbols), we obtained a quite good agreement with the val¬ 
ues of C computed from the mode eigenfunctions in Sect. 
14.1.21 at the exception of nearly pure-g modes for which C 
slightly overestimates (. As described in Appendix 14.21 we 
found that the agreement could be improved by releasing 
one of the approximations made by G13. We therefore pro¬ 
pose a slightly modified expression for ( given by Eq. IB.41 
(see Appendix [B] for more detail). The C computed with 
this latter expression are in excellent agreement with the 
actual C parameters (red symbols in Eig. [S]). This new ex¬ 
pression for C was tested on several other models of clump 
stars and systematically gave an agreement as good as the 
one shown in Fig. El We thus confirm that an estimate 
of the mode trapping can be obtained from the frequency 
modes without having to compute a seismic model of the 
star, as was claimed by G13. 

In principle, the rotational kernels of the modes are 
required to extract information about the rotation pro¬ 
file from the mode splittings. However, G13 showed that 
the rotational splittings can be related to the parameter C 
through a linear dependence whose coefficients depend only 
on the mean rotation rate in the g-mode cavity (Hg) and 
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Fig. 7. Rotational splittings as a function of the parameter which is an estimate of the mode trapping (see text). The circles 
correspond to the observed splittings (black symbols indicate modes that passed both frequentist and Bayesian significance tests 
and gray symbols modes that passed the frequentist test only). The black dot-dashed line shows a linear regression of the observed 
splittings as a function of and the gray shaded area indicates the uncertainties on the regression parameters. The dashed lines 
indicate theoretical splittings for 2-zone models that reproduce the observed rotation in the g-mode cavity and with flcore/S^surf 
of 1, 2, 10, or 100 from top to bottom (same colors as in Fig. HI). 
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Frequency (mHz) 


Fig. 6. Comparison between the mode trapping parameter 
of the 1 = 1 modes of the best-fit model for K1C7581399 com¬ 
puted from the mode eigenfunctions (Eq. 1151 black solid line) 
and the approximate expression ^ computed from the mode fre¬ 
quencies using the prescription of G13 modified in Appendix [B] 
(Eq. IB. 41 red filled circles and red dashed line). For reference, 
we also show the ^ computed with the original prescription of 
G13 (gray dotted line and open circles). 


the mean rotation rate in the p-mode cavity (fip): 

- (f}p)/(27r)) C + (f^p)/(27r) (17) 

Estimates of (fig) and (flp) can thus be obtained by per¬ 
forming a linear regression of the function Sv((). This was 
already applied to early red giants by D14, who found that 
it provided results in excellent agreement with the inver¬ 
sions performed from stellar models. However, they had 
been using the trapping parameters computed from stel¬ 
lar models. We here tested the method when replacing the 
parameters C by their approximate expression computed 
with Eq. IB.41 For KIC7581399, such an approach yielded 

(Hg)/(27r) = 164 ± 14 nHz 
(Hp)/(27r) = 88 ± 19 nHz 
(Hg)/(Hp) = 1.9 ±0.4 

These results are in very good agreement with the values 
obtained the mode eigenfunctions of the best-fit model in 
Sect. 14.11 This validates the simplified approach proposed 
by G13. 

4.3. Constraints on the core-envelope contrast for seven 
Kepler secondary clump star 

To estimate the amount of differential rotation for the seven 
targets selected in our sample, we have used the simplihed 
approach of G13 that was described and validated with the 
test-case of KIC7581399 in Sect. 14.21 For each star, the 
trapping _pf the mod e wa s estimated by computing the pa¬ 
rameter C using Eq. IB.41 Fig. 0 shows the mode splittings 
that passed the statistical test described in Sect. 12] as a 
function of As was already observed for KIC7581399 
in the previous section, it is striking to see how little the 
splittings vary with the mode trapping. For one star of the 


sample (KIC7467630), the splittings of p-dominated modes 
are found to be larger than the splittings of g-dominated 
modes. 

We then used Eq. [T7|to obtain estimates of (Hg) and 
(Hp) for all the stars of the sample. The results are given 
in Table 121 The ratios {Ctg)/{flp) were found to vary from 
1.3 to 3.2 for all the targets. This clearly conhrms that the 
radial differential rotation in secondary clump stars is much 
weaker than it is for RGB stars, as was suggested by the 
analysis of KIC7581399. However, we can exclude the sce¬ 
nario of a SB rotation for six of the seven selected targets, 
the last one (KIC7467630) being marginally consistent with 
a SB-rotating profile. 


5. Spectroscopic observations 

To complement seismic estimates of the rotation profile, 
we tried to obtain constraints on the surface rotation rate 
using high-resolution spectroscopic observations. For this 
purpose, we have observed three stars of the sample with 
the High Efficiency a nd Resolution Merca tor Echelle Spec- 
trograph(HERME^ Raski^,et_abL 201lL.Raskin & Van 


Winckel l2014l ) mounted on the 1.2 m Mercator Telescope 
at the Spanish Observatorio del Roque de los Muchachos 
of the Institute de Astrofisica de Canarias, with a resolv¬ 
ing power of Rhermes = A/AA~ 85000, where A and AA 
are the wavelength and the width per wavelength bin, re¬ 
spectively. The raw spectra were reduced and wavelength 
calibrated through Thorium-Argon reference spectra with 
the current version of the HER MES data reduction pipeline 
(Version 5. lR,askin et aLll^OllD . The radial velocities (RV, 
Table]?]) were obtained through a weighted cross-correlation 
of the wavelength range between 478 and 653 nm of each 
spectrum with a discrete G2 te mplate (jRaskin et al.l]2?)Tl] . 
iRaskin fc Van Winckell [2014] ). The post-processing and 
renormalization o f the spectra were performed following 
iBeck et ^ (|2015D . 

The total line broadening LB was determined from the 
equivalent wi dth of seven unbl ended Eel lines, following th e 
procedure of I Gray] (120051 ) and iHekker fc Melended (|2007D . 
The average broadening and its uncertainty, reported in Ta- 
ble0for the three stars were derived from the average and 
standard deviation of the individual total line broadening 
of the seven unblended Eel lines and deconvolved with the 
instrumental broadening of 1.76 kms“^. The stellar line 
broadening consists in contributions of the projected sur¬ 
face velocity, v sin i and the macro turbulence, Wmacro such 
that 


LBsteiiar sin i {vsmi)‘^ + V. 


2 

macro' 


(18) 


Disentangling the two contributions is challenging for red 
giant stars, as they have very narrow lines and Umacro 
was shown to have substantial contributions for this type 
of stars. Following_^he_scaling_relations_^on^Hekker & 
Melendez ( 2007! ) for the value of we estimated the 

contributions of both parameters to be on the order of 
~ 6 kms“^. However, owing to the large scatter in the 
data from which the scaling relations were obtained, we re¬ 
frain from deriving a usini for each star. Instead, we use 
the stellar line broadening as an upper limit to the value of 
usini by assuming that Wmacro = 0. 

In the present case, the inclination angles i of the stars 
can be estimated from seismology. For each star, the ex¬ 
traction of each rotational splitting provided an estimate 
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Table 7. Spectroscopic results for the three targets that were observed with the HERMES spectrograph 


KIG 

RV 

(kms“^) 

LBgteiiar sin i 
(kms“^) 

inclination i 

n 

LBgteiiar 
(kms“^) 

upper limit on Hsurf/( 27 r) 
(nHz) 

4659821 

-35.32 ±0.01 

8.5 ± 0.5 

79 ±4 

8.6 ± 0.6 

232 ± 25 

9346602 

-11.90T0.01 

8.4 ± 0.5 

83 ±2 

8.5 ±0.5 

158 ±21 

7581399 

-23.81 ±0.01 

8.9 ± 0.7 

79 ±5 

9.1 ±0.9 

193 ±29 


of i (see Tables |3] and n. By bringing these estimates 
together, we obtained estimates of i for the three stars, 
which are listed in Table [T] We used these values to cor¬ 
rect the stellar line broadening from the inclination effect, 
and thus obtained upper limits to the surface velocity of 
these stars. Since estimates of the radius of the three stars 
were already obtained from seismic scalings (see Tabled), 
we could translate these velocities into upper limits for the 
surface rotation rates, which are given in the last column 
of Table [71 These upper limits are of the order of the core 
rotation rates found through seismology (see Table 0, so 
the spectroscopic observations can be used only to rule out 
that the surface might be rotating faster than the core, but 
unfortunately not to test the core-envelope ratio that was 
obtained from seismology in Sect. H) 


6. Discussion and conclusion 

The aim of this paper was to measure the amount of differ¬ 
ential rotation in secondary clump stars (intermediate-mass 
core helium burning stars) from seismology. This required 
to measure the rotational splittings of both g-dominated 
and p-dominated modes from the oscillation spectra of 
these stars. This was a challenge because the splittings 
of clump stars are comparable to the mode widths, partic¬ 
ularly for p-dominated modes whose lifetimes are shorter. 
We followed both a frequentist approach using a maximum 
likelihood (MLE) method and a Bayesian approach with a 
Markov chain Monte-Carlo method to select only the sta¬ 
tistically reliable splittings. We found a sample of seven 
Kepler secondary clump stars with both g and p-dominated 
reliable splittings. 

To measure the amount of differential rotation in the se¬ 
lected stars using the mode splittings, we used an approach 
based on asymptotic theory proposed by G13, which has 
the advantage of being model-independent. We first vali¬ 
dated this method on secondary clump stars by confronting 
it to the results of classical inversion techniques for one tar¬ 
get of the sample (KIC7581399), and we proposed a slight 
modification to it. We applied this method to the selected 
stars and found evidence for a weak differential rotation in 
six out of the seven stars of the sample, with ratios between 
the core and envelope rotation rates ranging from 1.8 ± 0.3 
and 3.2 ± 1.0. The last target was found to be marginally 
consistent with a solid-body (SB) rotation profile. We tried 
to complement the seismic measurement of the internal ro¬ 
tation with spectroscopic measurements of the surface ve¬ 
locity. High-resolution spectra were obtained for three of 
the selected targets, from which we deduced upper limits to 
the surface rotation rate. These limits confirmed that the 
surface cannot spin faster than the core, but they were not 
sufficient to test the amount of differential rotation found 
with seismology. 


Our results clearly show that secondary clump stars 
have a much milder differential rotation than their RGB 
counterparts, for which core-envelope rotation ratios above 
20 were observed (e.g. G13) . Following the approach of 
iTavar fc Pinsonneaultl (120131 1. we compared the total AM 
of the secondary clump stars of our sample to the estimated 
AM at the terminal-age main sequence (TAMS). Assuming 
SB rotation at the TAMS and typical surface rotation rates 
(50 to 150 kms“^), we expecljj a total AM at the end of 
the MS around (4 ± 2) 10"^® erg.s. By comparison, using 
the mean core and envelope rotation rates derived in this 
study and our optimal model of KIC7581399 as a reference 
model for secondary clump stars, we obtained estimates of 
the total AM on the clump around 2 10"^® erg.s for the seven 
stars of the sample. This value is of the same order of mag¬ 
nitude as the total AM expected at the TAMS, and the 
slightly lower AM found on the clump might be caused by 
mass loss, which is expected to occur between the TAMS 
and the triggering of core-He burnin g (not taken into ac¬ 
count bv lTavar fc Pinsonneaultl^2013^ . 


The weak radial differential rotation that we find for 
secondary clump stars means that a very efficient redistri¬ 
bution of AM occurs either before or right after these stars 
settle on the clump. If this strong coupling occurs before 
He-burning is triggered in the core, it has to operate on 
a very short timescale since these stars cross the subgiant 
phase on a thermal timescale. To illustrate this point, a 3- 
Mq star evolves from the MS turnoff to the clump in only 
15 Myr, which amounts to about 5% of the time it spends 
on the MS. For comparison, it takes around 700 Myr for a 
I.S-Mq star, i.e. about 34% of the time spent on the MS. 
So if the same mechanism is responsible for AM redistribu¬ 
tion in low-mass and intermediate-mass stars, the present 
results provide more stringent constraints on its timescale. 
The other possible explanation is that AM could be redis¬ 
tributed after the triggering of He-core burning. Internal 
gravity waves (IGW) excited both at the bottom of the 
convective envelope and at the top of the convective core 
generated by He- burning could efficiently couple the core 
to th e envelope llTalon fc Gharbonnell 120081 . [Fuller et al.l 
I 2 OI 4 I I. iMaeder fc Mevnetl ( 2014| l have also recently shown 
from order-of-magnitude estimates that a fossil magnetic 
field attached to the core could produce efficient coupling 
between the core and the envelope during the He-burning 
phase. To determine at which stage the coupling occurs 
would require to find intermediate-mass stars in the sub¬ 
giant phase or at the bottom of the RGB. Even though this 
phase is very fast, there might be a handful of such targets 
among the Kepler targets. 


® To derive this estimate, we used the moment of inertia of 
the optimal model derived for KIC7581399 (Sect. 14.1. HI whose 
evolution was stopped at the TAMS. 
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Regardless of when the AM redistribution occurs, it has 
to happen on a very short timescale, which should provide 
valuable observational constraints to theoretical models of 
AM transport mechanisms in intermediate-mass stars. It is 
important to note that even though an efficient AM trans¬ 
port is needed, it does not generally lead to SB rotation 
since the SB rotation case was significantly ruled out for six 
of the seven targets of our sample. Interestingly, the two 
targets that were found to have the lowest core-envelope ro¬ 
tation contrast are the ones with the largest values of AIIi. 
Sinc e AIIi increases dur ing the core He-burning phas^l (see 
e.g. iMosser et al.ll20l3 l. this could mean that SB rotation 
progressively builds up during this period. More data are 
required to confirm this hypothesis. In this con text, the ob¬ 
serva tions of the selected space mission Plato (iRauer et al.l 
120141 1 will be particularly useful. 
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Appendix A: Monte Carlo 

To verify that AA is distributed as a with 2 degrees of 
freedom under the Hq hypothesis, we performed a Monte 
Carlo simulation. We generated a mode profile with fre¬ 
quency, height, and linewidth typical of oscillation modes 
of stars in the clump. For each iteration, we included noise 
following the distribution of a with 2 degrees of freedom 
to simulate an observed power spectrum. Two fits were 
performed, either under the Hq hypothesis (best likelihood 
£ 0 } or under the Hi hypothesis (best likelihood £ 1 ), and 
we stored the value of AA = —2 [In(^l) — In(^o)] for each 
iteration. Fig. lA.ll shows the distribution that we obtained 
for AA. For comparison, the distribution of a with 2 de¬ 
grees of freedom is overplotted (red). We confirm that the 
two distributio ns are i ndeed very similar, which validates 
the assertion of lWilksI (llQSST l in our particular case. 
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Fig. A.l. Distribution of AA obtained with a Monte Carlo 
simulation with 500 iterations (black curve). The distribution 
of a with 2 degrees of freedom is overplotted (red curve). 


Ap pendix B: Modification of the lGoupil et al. 
( |2013 ) formula 


obtain estimates of Au, Atti and q in Sect. we u sed 
Eq. d] and 0 for 9p and 9g respectively in Eq. IB.21 in 
order to be consistent. 

— Finally, G13 assumed that the factor exp(0e) in Eq. IB.21 
was close to 1 owing to the narrow evanescent zones in 
red giants. As mentioned above, this factor corresponds 
to (4q)“^/^, where q has been measured from the mode 
frequencies. We found that this factor is indeed not far 
from unity (it varies from 0.7 to 1.1 for the stars of the 
sample) but taking it into account slightly improves the 
agreement with the actual values of ((. 


In the end, we obtained the following modified expres¬ 
sion for ^ 


c 



COS^ 


1 - 'll 

1/2 ABi 

i/AHi ^Sj 

COS^ 

> 1 

V 

< 


-1 


(B.4) 


This expression was used to estimate the mode trapping in 
a model of KIC7581399 in Sect. I4.2l and an excellent agree¬ 
ment has been obtained with the C parameter computed 
with the mode eigenfunctions and Eq. [T3(see Fig. O. 


We start by briefly recalling the method proposed by G13. 
By using approximate expressions of the eigenfunctions in 
the p-mode and g-mode cavities through asymptotic anal¬ 
ysis, they approximated the mode trapping parameter ^ by 
the quantity 


c 





-1 


(B.l) 


where dp and dg were defined in Eq. [31 The coefficients 
a and c are the amplitudes of the eigenfunctions in the g- 
mode and p-mode cavities, respectively. The matching of 
the solutions in the evanesce nt zone provides th e following 
relation between a and c (see lUnno et al.lll989ri : 


c 

a 


o COs(dg) 

cos(dp) 


exp (de) 


where 

r° 

de = / ndr 
J i-b 


(B.2) 

(B.3) 


and in the evanescent zone. We note that the 

factor 2exp(de) corresponds to where the parameter 

q was introduced in Eq. [3] and whose value was already 
estimated from the mode frequencies for all the stars of the 
sample (see Sect. 12 .21) . We here made minor modifications 
to the expression of C proposed by G13 by releasing three 
approximations that they made. 

— First, in their equivalent of Eq. EJ G13 argued that 
cos dp Ri 1 since dp k, Upir. Technically, this is true only 
when the modes are p-dominated and should therefore 
not be assumed for all modes. We kept this term in the 
expression of C- 

— Secondly, the approximate expressions for dp and dg 
used by G13 (Eq. A12 and A18 of their paper) in the 
phases of the cosine functions in Eq. IB.21 neglect the 
phase shifts £p and £g. Since we have included them to 
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